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Abstract. We consider radiation transport theory applied to non-dispersive but 
refractive media. This setting is used to discuss Minkowski’s and Abraham’s 
electromagnetic momentum, and to derive conservation equations independent of the 
choice of momentum definition. Using general relativistic kinetic theory, we derive 
and discuss a radiation gas energy-momentum conservation equation valid in arbitrary 
curved spacetime with diffractive media. 


PACS numbers: 41.20.Jb, 05.60.-k, 42.15.-i 

1. Introduction 

Radiative transfer is a mature area of research, with important application both in 
laboratory and astrophysical systems (see, e.g., [1,2] and references therein). Much 
is thus known about the properties of the relevant radiation transport equations in 
non-dispersive, non-diffractive media in flat spacetime. However, less is known in the 
case of dispersive and/or diffractive media, and the addition of spacetime curvature 
reduces the number of relevant publications further (for a representative but incomplete 
selection see [2-8]). Different ways to tackle the problems involved in formulating a 
general theory of curved spacetime radiation transport in diffractive media exists in the 
literature. In many cases, the treatment of moments of the transfer equations, and their 
respective conservation equation, is lacking, as is the discussion of the interpretation 
of the radiation fluid moments. The latter can be viewed as a nontrivial task, since a 
clear-cut definition of the radiation momentum density is still under discussion (although 
claimed otherwise by some authors, see, e.g., [9]). There are exceptions in the literature 
though. For example, Anderson & Spiegel introduces an optical geometry on top 
of the curved spacetime, by introducing an effective metric which incorporates the 
refractive index of the medium (similar to Gordon [10]). In this effective geometry, the 
transfer equation takes the standard form, and a treatment of the fluid moments is thus 
straightforward, although their interpretation is somewhat obscured by the presence 
of the effective geometry. Moreover, the proper conservation equations for the fluid 
moments is not derived. 

As noted by many authors, related to the problem of radiative transfer in 
diffractive media is concept of photon momentum and the Minkowski-Abraham debate. 
Although the definition of Abraham is preferred due to a number of theoretical reasons 
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(symmetries [9, 11], derivations from microscopic theory [12] etc., see [13] for an 
overview), measurements have not given a dehnite answer [14,15] (for a review, see [16] 
and references therein). Due to the problem of separating the electromagnetic from 
the material degrees of freedom [17,18], the problem of defining the electromagnetic 
momentum in a refractive medium has persisted. However, it seems reasonable that 
the Minkowski momentum should be treated as a pseudo-momentum, partly depending 
on material contributions, while the Abraham momentum is a proper electromagnetic 
momentum (for a discussion, see [19]). 

Here we will analyse the radiative transfer equations and derive macroscopic 
conservation equations in terms of both Abraham’s and Minkowski’s definitions. It 
will also be shown that variables can be used as to remove the problem of momentum 
definition when deriving conservation equations. Furthermore, a completely general 
energy-momentum conservation equation for a radiation fluid in curved spacetime with 
a refractive medium is derived and discussed. It is shown that the combination of 
curvature and refraction gives important contributions to the conservation equation, 
even if the refractive index is spacetime homogeneous. The equation is discussed in the 
context of cosmological models. 


2. Ray dynamics and kinetic theory 


In mechanics, given the action S' of a system, we may define the momentum and 
Hamiltonian according to 


Pa = 


dS 

dx^ 


H- 


( 1 ) 


respectively. From this, the equations of motion for a single particle can be written in 
terms of Hamilton’s equations, i.e. 


X 


a 


Pa 


dH 

dpa 

dH 

dx'^ 


(2a) 

m 


In geometric optics [11], we introduce the eikonal 0 = 0(t,a;“) by writing the 
electromagnetic field F in the form 

F = aexp(i0). (3) 


From the eikonal, we define the the wave-vector and frequency of the field as 


kr, = 


d(j) 


UJ = -- 


d(j) 


dx’^ dt ’ 

in analogy with (^. The equations of ray optics becomes 
duj 


x^ = 


dkn ’ 


k 


a 


duj 
dx^ ’ 


( 4 ) 

(5a) 

(56) 
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and from this also follows that duj/dt = duj/dt. These equations are the mass-less 
quasi-particle analogous of ©• It is also clear from the eikonal dehnitions that the 
natural variables describing the photons are the frequency and the wave-number rather 
than the energy and the momentum. 

Photons moving in a isotropic non-dispersive but refractive medium will satisfy a 
dispersion relation of the form 
ck 

uj = — 
n 

where, in general, the refractive index n is spacetime inhomogeneous. The equations of 
motion m for individual photons in such a medium then becomes 

(7a) 


k 


a 


n 

oj dn 
n dx'^ 


(76) 


Before proceeding further, a discussion of the dehnition of the photon momentum 
is in place. We note that we have chosen as our fundamental variables and ov, k^- 
The relation to Hamiltonian particle dynamics becomes obvious if we choose H = hx 
and Pa = hka, and it is therefore tempting to assume this direct relationship. However, 
in electromagnetic theory, there are two distinct ways to define the momentum of the 
electromagnetic field. According to Minkowski, the momentum is proportional to D x B, 
while according to Abraham it is proportional to E x H (see Ref. [9] for a discussion). 
Furthermore, in measurement of photon momenta in isotropic diffractive media, there 
are two distinct forms as well, namely p = hk = hwn/c and p = tux/nc = hk/n^. 
We note that the former is in agreement with Minkowski’s momentum density and a 
number of experiments, while the latter dehnition is consistent with Abraham’s choice 
in terms of the Poynting hux (see [16] and references therein), which also follows from 
detailed microscopic considerations as well as symmetry arguments. On the other hand, 
measurements indicates that p oe. n, consistent with Minkowski’s momentum dehnition 
in a dielectric medium. The discrepancy between the two dehnitions can be attributed to 
the contributions of the medium in Minkowski’s dehnition, while Abraham’s dehnition 
(used here) is the ‘proper’ momentum of the photon [19]. However, we note that while 
the form (IZ3) of the group velocity is valid independently of the momentum dehnition 
(as expected), (17^ formulated in terms of the momentum takes diherent forms for the 
Minkowski and the Abraham dehnition. 

Given a spectral distribution Ji^{t,x°‘,ka) of photons, the absence of collisions 
dehnes the on-shell Vlasov equation for t/F according to 


^ + ka^^ = 0 , 


( 8 ) 


dt dka 

expressing the phase space conservation of quasi-particles. 

Using the spectral distribution function dehned on-shell, i.e. with satished, 
we can now dehne the macroscopic observables, from which a huid theory can be 
constructed. 
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Next, we define macroscopic variables as moments of the distribntion fnnction In 
general, we define ('0) = (/ ^dfk)~^ f xf^dfk to be the statistical average of the 
fnnction f), which may be dehned over the fnll phase space. Moreover, in these dehnitions 
we closely follow Ref. [1]. 

— The number density N is dehned according to 


N = j ^d?k. 

— The energy density fi is dehned as 

R — J kw^d^k. 

— The average fluid velocity is given by 
and from this we may also dehne the thermal (or random) fluid velocity tc" nsing 




(9) 

( 10 ) 

( 11 ) 

( 12 ) 


snch that (w“) = 0. 

— The radiation energy flux q°‘ is dehned snch that g" dSa is the rate of energy 
how across the snrface dS^. Thns 


q°‘= hujw'^^d^k. 


(13) 


— The momentum density IT" of the photon in a refractive medinm can be 
formnlated in two ways, in accordance with Minkowski or Abraham, and they read 

,2 


= / hk^jRd^k = / ^hwx^jRdt^k, 


and 


= 


hk"^ 




■^d^k = / —hujx'^^d^k, 


(14a) 


(146) 


respectively, where we have made nse of ®-ca. 

— The pressure tensor P“^ can similarly be dehned in two distinct ways. Since the 
pressnre tensor is dehned to be, for an observer comoving with the average hnid 
how, the rate of transport of the a component of momentnm per nnit area of a snrface 
orthogonal to the [3 coordinate (or tetrad) axis, we have 

,2 


p^ = / hk^w^^d^k = / ^hxw^w^^d^k, 


and 


p«/3 _ 




-hk^w^^d^k = / -^hujw^w^^d^k, 


(15a) 


(156) 


for the Minkowski and the Abraham momentnm, respectively. Here we have again nsed 
ds-cg. and we note that the pressnre tensor is symmetric. 
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Furthermore, the energy flux and the momentum density is related to each other 
through the relations 
nf 

+ g“), (16a) 

and 

n^ = 4(h«“ + g“), (166) 

where we have used general relation + w°')/(?. 

Higher order moments can be deflned analogously, but the above will suffice for 
our purposes, and with these definitions we are ready to set up the fluid equations. We 
note that the macroscopic quantities defined above are valid also for dispersive media, 
i.e. when the refractive index depend on the frequency uj or, in the anisotropic case, the 
wave vector 


4. Fluid equations 


A hierarchy of fluid equations may be obtained from (jH]) by taking the moments with 
respect to suitable microscopic quantities. Given a microscopic variable (which may 
be a tensorial object), the general fluid conservation equation takes the form 

+ A (Ar(±“^)) = (17) 

where is defined in accordance with (jH)). Thus, Nipp) and N{x°''ijj) represent the 
macroscopic density and macroscopic current, respectively, of the microscopic variable 
■0, while the right hand side of (ini) acts as a source/sink. Finally, the fluid hierarchy of 
equation can be terminated by assuming a set of thermodynamic relationships between 
the macroscopic variables, such as an equation of state. 

Putting = 1, we obtain the conservation equation for the number of quasi¬ 
particles 


ON d 

"aT ^ ^ “ ' 


(18) 


The energy conservation equation is obtained by taking tp = hx. Thus, for an 
observer comoving with the average radiation fluid flow we obtain 


A + („„. + „.) = _#^ hi. 

dt dx^ n dt 


(19) 


Similarly, the momentum conservation equation is obtained by '^ = and takes 
the form 


an 


dt 


Ma ^ d 


dxf^ \ 


I 


+ !!:„ + P« 1 - 


n dx° 


(20a) 


or 

dt dxd 



fi dn 
dx°‘ 


n 


Acx. 


n 


d . d 


dxd 


n 


n \c 


- + P 

m \ 


Aa 


dn 
dxd ’ 


{20b) 
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depending on if we use (ina) and (USD or (I146jl and ()156jl . respectively. We furthermore 
note that ()2()a|l and can be written in a slightly more symmetric form, using the 
relations according to 


an' 


1% a /n- „ 

dt 




+ 


2n^ 




I pa/3 


= 

n dxd 


dn 


(21a) 


and 

dt 


a 

dxd 



= ^ ya/3 dn _ 2 dn 
dxd n ^ dt 


2 


n 





dn 

dxd 


( 216 ) 


respectively. It is straighforward to see that (l2Ual) - (l206|) can be obtained from each 
other using the relation = n^TVf. 

The system of equations presented above can be closed if we choose a 
thermodynamic relationship between certain quantities. For an close-to-equilibrium 
system, the pressure tensor becomes nearly isotropic, and we can write ~ 
where P = From (jlbojl and (1155jl . we obtain 

Pm = |/i Pa = (22) 


respectively. 

The remaining freedom in the equations is removed by choosing the observer. We 
will here use the particle frame, where = 0. 

Using these choices, we obtain a set of equations independent of the choice of 
momentum dehnition, in terms of the energy density and heat flux, according to 

dt dx^ ndt' ^ ^ 

and 


dt 3n^ dxd n^ dx^^ n dt ’ 

from (HHii)- (dH), (unD, and (Ed-iEd, respectively. We may from this derive general 
wave equations for the energy density and energy flux for radiation in diffractive media. 


5. Covariant conservation laws 

In this section we will make use of the 1 + 3 orthonormal frame (ONF) approach 
(see [20] for an overview), since it allows for a general relativistic treatment, simplihes 
calculations, and gives the natural Cartesian like reference frame for an observer moving 
with the timelike frame direction. All vector and tensor quantities will be projected onto 
this frame. 

We dehne a set of frame vectors {e^}, a = 0,...,3 such that Sa ■ ei, = rjab gives 
the constant coefficients of the metric: Pab = Pab = diag (—1,1,1,1), i.e., a Lorentz 
frame. The Ricci rotation coefficients T%c, antisymmetric in the hrst two indices, giving 
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the kinematics of spacetime are defined by 'V^ei, = T%cea. Here Vc is the covariant 
derivative with respect to the Lorentz frame. 

The distribution function ^ will now be a function of the canonical phase space 
variables and ka, and the conservation of phase space density can be written 


dJR 

~dX 


L[A 


+ ka 


dJ" 

dka 




(25) 


where the over dot stands for d/dX, where A is an affine parameter along the photon 
path, and we have introduced the Liouville operator L = x^'Sa + kad/dka- Furthermore, 
for the sake of generality, we have added the collisional operator representing photon 
emission and absorption If /“ are the external forces acting on the pencil of light, 
the covariant derivative along the photon path is 
Dku 

-^ = h- nakcX^ = fb (26) 

dX 


With the Hamiltonian H = H{x°',ka) in eight dimensional phase space, the equations 
of motion thus reads 

dH 


x^ = 


dka 


ka = -ea{H) + r 


(27a) 

{27b) 


generalising ()5aj) - (l^ . Here the last term of ()27^jl gives the gravitational contribution 
to the equations of motion. We denote the observer four-velocity by ?7“, normalised 
such that IR^Ua = —c^. We partially £x the frame by letting the observer four-velocity 
U°' = 6q, and split spacetime quantities with respect to U°‘. The space metric orthogonal 
to takes the form 


hab = Vab + UaUb/c^, (28) 

and the wavevector can be written 

ka = OvUa/c^ + kia, (29) 

where u = —U°‘ka, k = {h^'^kakf}^/"^ , and C = hfkiy/k. We note that i^-ia = 1 and 
£^7/“ = 0. In these variables, the dispersion relation © retains its form. 

There is a certain arbitrariness in the choice of Hamiltonian. We choose the 
Hamiltonian H such that H = 0 gives the dispersion relation (and OH/dka 7^ 0 on 
the dispersion surface), and the variables x“,a;, /c, and ^a is treated as independent, with 
the vanishing Hamiltonian as a constraint. Suppose the dispersion relation takes the 
form 


UJ = W{x^,k,la)- (30) 

f The process of emission and absorption is thoroughly treated in any book on radiation transport, 
see e.g. [1], and will therefore not be discussed further here, apart from the following brief comment. 
As an example, in case of applications to the early universe, the major contribution to the collisional 
term would be in the form of Thomson scattering, but there are of course numerous other scattering 
events that could dominate in other parameter regimes. 
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H = -uj + W{xfk,ia) 

satisfies the necessary criteria listed above, giving the equations of motion 


( 31 ) 


i:“ = + 


dW 


via 

and 


dka 

ka = -ea{W) + 

, and, moreover, H = eo{H) = eo{W). With (P), W{x°',k) = kc/n{x°‘), 


(32a) 

(326) 


dW 

dka 


ct 

n 


00 


eaiW) = - ea{n). 


n 


(зза) 

(ззб) 


Thus, with the particular dispersion relation (P, the kinetic equation (ES) becomes 


c£“ 


n 


00 


H-) ea(^) + —Gain) + T^abkcX 


in 


J dka 


= ^. 


In the literature, the common definition of the energy-momentum tensor is 


lab 


=fi / K, 


(34) 


(35) 


where K. — Idet^l /bj) d'^k is the invariant volume measure in momentum space. 

This definition is consistent with the Minkowski definition of the photon momentum in 
diffractive media. Furthermore, in non-diffractive media, this is a conserved quantity, 
due to the one-to-one correspondence between and ka- However, in diffractive media, 
this correspondence is lost, since 


i;“ = c2— + -(l-n2)r= - 
00 n 00 


9 


ab 




h 


.^ab 


kh 


(36) 


[see and (I33a|) ]. Thus, in diffractive media, the energy-momentum tensor may 

alternatively be defined according to 

h 


Tf = h 


2 ’Cl * /’ 

00 X X 


o/F K = 


oox^^x^jR 6{H) d^k, 


(37) 


J C^ c c c^ 

in accordance with Abraham’s photon momentum definition. We note that when n —>■ 1, 
we regain the common definition of the energy-momentum tensor in non-diffractive 
media. Moreover, since (jsni) holds in general, we have the relation 


rpab _ 

■^A — 


9^ 


4 1 


9d 


4 ] h’’ 




T; 


cd 

M 


(38) 


between the two definitions of the energy-momentum definition. 

From any energy-momentum tensor, we may define the fluid quantities used in the 
preceding sections. We thus have the relativistic energy density fi = TabU°'U’^ with 
respect to f/“, the relativistic momentum density n“ = the scalar pressure 

P = (c^Taft/i“^)/3, and the full anisotropic pressure Pab = Tcdh^ah%- By construction n“ 
and Pab are spacelike quantities, i.e., orthogonal to U°'. 
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The equations of conservation of energy and momentum can be obtained by 
multiplying (El by /c^, and integrating over the momentum space variables on 

the dispersion surface. The result is the energy-momentum conservation equation 


v„rf = 2 


1 


ad 


n^g' 


-n^(l 


— 1 


-id)e 


..ab 


9 


4 > ^ 


ab 


UdUe + 2U^Udg\ - 4h(' 




-t- C\ 
n 


(39) 


where = f (huj^x^/c‘^)'^ k is the collisional contribution, and Va is the covariant 
derivative with respect to the Lorentz frame. We note that as n —>■ 1, all terms 
on the right hand side, except C^, vanishes and we obtain the standard conservation 
equation for non-diffractive media. Moreover, using the transformation we obtain 
the corresponding equation in terms of the Minkowski variables. The hrst term on 
the right hand side gives the coupling between spacetime kinematics and the media 
refractive properties. Thus, when we have either flat spacetime or non-diffractive media, 
this term vanishes. The second term arises from the spatial and/or time-like dependence 
of the refractive index. Equation ()d9|l gives the evolution of energy and momentum of a 
radiation fluid in diffractive media on an arbitrary curved spacetime. Thus, El is the 
proper starting point for analysing photon gas dynamics in general relativistic gravity. 
As compared to the case of a Minkowski background spacetime, there is a signihcant 
alteration to the equation due to the novel coupling between gravity and the diffractive 
media, as given by the hrst term on the right hand side of Ei- 

As a simple example of curvature effects, we use El to derive the appropriate 
energy conservation equation in a Friedmann-Lemaitre-Robertson-Walker spacetime, 
i.e. spatial homogeneity and isotropy is assumed. The spacetime is characterised by the 
time dependent scalar quantities fv (energy density), P (pressure), 0 (expansion), and 
n (refractive index) (see [20] for details). Thus, El gives 


dpi 4 3pidn 

dt 3 ~ ^ n dt' 


(40) 


where we have used Pa = p./'in^ [see (El], and d/dt = Of course, the analysis 

of the equation requires a specihed density dependence of the refractive index, and 
therefore also a modihed equation of state (for an example, see [21,22]), and will not be 
pursued further here. 

We note the possibility to formulate El in terms of an effective energy-momentum 
tensor in certain cases, in particular for homogeneous spacetimes [23]. For example, we 
may dehne pes = 'nfp, so that dH takes the standard form for this new effective density. 
In general though, this is not a consistent approach [23], and can apparently only be 
overcome by introducing an effective geometry (see, e.g., [7]). 


6. Conclusions 

We have discussed the consequences of the different electromagnetic momentum 
dehnitions, due to Minkowski and Abraham respectively, in the context of radiation 
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fluid dynamics. Starting from a kinetic description, a set of flnid eqnations was derived 
and compared for the different definitions of flnid variables. It was fonnd that by 
expressing the equations in terms of certain variables, they became independent of the 
choice of momentnm dehnition. Finally, from a general relativistic kinetic theory the 
energy-momentum conservation equation valid for a radiation gas in a refractive medinm 
on an arbitrary spacetime, including collisional effects, was derived and discussed in a 
cosmological setting. 
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